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Abstract 

i , 

jj2 \ Argon with an admixture of CF 4 is expected to be a good candidate for the gas mixture to be used 
"~J ' for a time projection chamber (TPC) in the future linear collider experiment because of its small 
q ■ transverse diffusion of drift electrons especially under a strong magnetic field. In order to confirm 
"5a ■ the superiority of this gas mixture over conventional TPC gases we carried out cosmic ray tests using 
a GEM-based TPC operated mostly in Ar-CF 4 -isobutane mixtures under - 1 T axial magnetic 
o , fields. The measured gas properties such as gas gain and transverse diffusion constant as well as the 
' observed spatial resolution are presented. 
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q 1. Introduction 



PACS: 29.40. Cs, 29.40.Gx 



A strong candidate for the central tracker of the future linear collider (LC) experiment [H, [i| is 
\ a large time projection chamber (TPC) [3]. The LCTPC is expected to have unprecedentedly high 
k> \ momentum resolution under a 3 - 4 T magnetic field and good two-track resolving power to precisely 
5_i ■ reconstruct high momentum lepton tracks and each of the charged tracks in densely packed jets. In 
^ order to fulfill these requirements a micro-pattern gas detector (MPGD) is planned to be used for 
the endplate readout device since a conventional multi-wire proportional chamber (MWPC) readout 
suffers severely from the E x B effect in the vicinity of its sense wire planes under a high magnetic 
field jij. A GEM-based TPC prototype has been proved to work stably in conventional gas mixtures 
and to give satisfactory spatial resolution (see, for example, Ref. j^j). 

However, the target spatial resolution of better than 100 fim in the r-cf) plane for high momentum 
tracks is still ambitious and difficult to achieve with conventional gas mixtures such as TDR gas (Ar- 
CH 4 (5%)-C02(2%)) or P5 gas (Ar-CH 4 (5%)), in particular at long drift distances (> 2 m) because 
of unavoidable diffusion of drift electrons even under a strong axial magnetic field. 
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Figure 1: Cross sections of argon, CH4, CF4 and CO2 as function of the electron energy See Refs. 13, 11 1 for 
each of the components contributing to the cross sections shown in curves. 



The spatial resolution in the pad row direction (<7x) as a function of the drift distance (z\ 
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where cxxo is a constant term, D is the transverse diffusion constants, and iV e fj the effective number 
of electrons^]. It should be noted that Eq. (1) contains no angular terms (the angular wire and the 
E x B effects) except for the angular pad effect, which is implicitly included in the constant term 
(oxo). The value of crxo becomes large as the track angle increases with respect to the pad- row 
normal. 

From our experience oxo is about 50 jum (or better) for tracks perpendicular to the pad row 
direction if the charge spread in the amplification gap is not too small compared to the readout 
pad pitch and the S/N ratio of the readout electronics is good enough. The parameter N e g has 
been measured to be around 20 in argon-based gases and for a pad height of ~ 6 mm 0, 0] , which 
is consistent with a simple estimate taking into account the primary ionization statistics and the 
avalanche fluctuation in the amplification gap [zj]. 

If we require the azimuthal resolution of 100 /mi at z = 200 cm the diffusion constant (D), which 
is essentially the only free (controllable) parameter depending on the choice of gas mixture, needs 
to be smaller than 30 /iin/^cm. The diffusion constant of drift electrons under the influence of an 
axial magnetic field (B) is given by 



D(B) = D(B = 0)/J1 + (lot) 



(2) 



1 The diffusion constant D is related to the diffusion coefficient (D*) through D 2 — 2D* /W, where W is the 
electron drift velocity. 

2 Equation (1) may be inappropriate for small z as will be seen in Section 4.2. Even in that case, however, it 
describes the asymptotic behavior of the spatial resolution at long drift distances. 
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where u = e • B/m, the electron cyclotron frequency, and r is the mean free time of drift electrons 
between collisions with gas molecules. Therefore we need a gas mixture in which D(B = 0) is small 
(cool) and r is fairly large (fast) under a moderate drift field (E)$ 

A guideline to obtain a cool and fast gas is suggested by the deep and broad Ramsauer minimum 
in the electron cross section in argon, which is located around the electron energy (e) of 0.2 eV. It 
is necessary somehow to confine drift electrons near the Ramsauer minimum in argon, where the 
electrons are relatively cool and the cross section is small (r is large), under a moderate drift field. 
Actually, argon itself is a very slow and hot gas and the drift electrons get hot quickly even under 
a weak drift field. Therefore it is necessary to add an efficient moderator, i.e. a small amount 
of molecular gas with its Ramsauer minimum at around e = 0.2 eV and with large inelastic cross 
sections at e > 0.2 eV in order to efficiently cool down relatively hot electrons jsl, [if. 

Good candidates for the additive gas are CH 4 and CF 4 . The cross sections of relevant gases to 
electrons are seen in Fig. [1] as function of the electron energy. Figure |2] shows the mean electron 
energy as a function of the electric field for some pure gases and gas mixtures given by Magboltz (lll|R. 
All the results of Magboltz (version 8.5) presented in the present work are for gases at NTP (20°C, 1 
atm). Examples of the electron energy distribution are shown in Fig. [3] for some gases of interest at 
several electric fields. The figure demonstrates that CH 4 or CF 4 molecules in argon are very efficient 
to keep the electron temperatures around those corresponding to the Ramsauer minimum, even under 
relatively high electric fields. Furthermore, CF 4 is superior to CH 4 as a moderator additive because 
of its much larger inelastic cross section to electrons with energies above the Ramsauer dip, due to 
vibrational excitation (see Fig. 1). 




Figure 2: Mean electron energy as a function of the electric field given by Magboltz for (1) Ar (100%), (2) Ar-CH4(5%), 
(3) Ar-CF4(3%)-isobutane(2%), and (4) CO2 (100%). Electrons are almost thermal in pure CO2 in the range of electric 
field shown in the figure. 

Although they are fast and (relatively) cool themselves it is favorable to dilute them with argon 
in order to lower the drift field corresponding to the drift-velocity maximum by effectively reducing 



3 The electron drift velocity is given by W = e ■ E/m ■ t with e (m) being the electron charge (mass). A large value 
of t, therefore, means a fast gas. 

4 The mean electron energy can be estimated experimentally from the characteristic energy (ejj): (e) ~ §£k = §■ , 
where D* is the diffusion coefficient and [i the mobility of electrons 12]. 
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Figure 3: Electron energy distribution in some gases given by Magboltz. The electric field of 80 V/cm (250 V/cm) 
approximately corresponds to the drift-velocity maximum for Ar-CH4(5%) (Ar-CF4(3%)-isobutane(2%)) while the elec- 
tric field of 35 V/cm (130 V/cm) corresponds to the diffusion minimum for Ar-CrLi(5%) (Ar-CF4(3%)-isobutane(2%)). 



their inelastic cross sections, and at the same time, the operating high voltage of MPGDs. A low 
drift field is favorable especially for a large TPC in order to limit the high voltage for the central 
membrane to a moderate level. It is worth mentioning here that the drift field should be chosen so 
that the drift velocity be as close as possible to its maximum since the gas pressure in the LCTPC 
follows atmospheric pressure. 

We have already tested several candidate gases, including TDR gas and P5 gas, for the LCTPC 
using a small prototype equipped with an MPGD (MicroMEGAS (l3| or GEM [lij]) in beam tests [j§, 

The resultant transverse diffusion constants were always close to those given by Magboltz under 
- 1 T magnetic fields. According to Magboltz the diffusion constant in P5 gas (TDR gas) at B 
= 4 T, and E ~ 80 V/cm (~ 220 V/cm) corresponding to the drift- velocity maximum, is about 40 
/jm/^cm (60 /im/^cm) [l5[, which does not meet our requirement mentioned abovej^ 

As for Ar-CF4, which is cooler and faster than P5 gas, we were not sure about the reliability 
of the Magboltz prediction although some measurements using prototype TPCs mainly equipped 



with MicroMEGAS exist [1614221] . We therefore conducted a series of cosmic ray tests using a small 
prototype TPC with a GEM readout operated in this possibly promising gas mixture of Ar-CF 4 (- 
isobutane) under a magnetic field of up to 1 T. In addition to the measurement of diffusion constant 
it was important as well to demonstrate successful operation of a GEM-equipped TPC in gases 
containing CF 4 since CF 4 molecules dissociatively capture relatively hot electrons in the vicinity of 
GEMs§ It should be noted that the loss of drift electrons at the entrance to the first GEM could 



5 A cool and slow gas such as DME (dimethyl ether) or CO2 (as in the case of TDR gas) is not a good additive 
since it makes the mean free time of the drift electrons (r) small because of a large cross section to low energy electrons 
(see Fig. [T]for the cross section of CO2). It certainly reduces the diffusion constant at B = T but, at the same time, 
makes the reduction of D under a magnetic field less effective. 

6 Significant deterioration of the energy resolution for soft X-rays due to the electron attachment in gases containing 



CF4 has been observed with proportional tubes 8, 23, 24] and with an MWPC-equipped TPC [25 
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be fatal to the performance of a TPC, i.e. the spatial and the dE/dx resolutions, since this loss can 
not be compensated by the subsequent gas multiplication. 

The experimental setup is briefly described in Section 2. The results of gas gain measurements 
are given in Section 3, while the measured drift properties of the gas and the performance of the 
prototype are presented in Section 4. Section 5 is devoted to the conclusion. 

2. Experimental setup 

A photograph of the prototype (MP-TPC) is shown in Fig. HI It consists of a field cage and 
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Figure 4: Photograph of the prototype just before installation into the gas vessel. 

an easily replaceable gas amplification device attached to one end of the field cage. Gas amplified 
electrons are detected by a pad plane at ground potential placed right behind the amplification 
device. A drift electrode is attached to the other end of the field cage. The maximum drift length is 
257 mm. 

In the present test a triple GEM is used as the amplification device. The triple GEM, CERN 



standard [26J, has 1.5 mm transfer gaps and a 1.5 mm induction gap. The high voltages applied 
across each of the GEM foils and the transfer and induction gaps are set equal with a resistor-chain 
voltage divider. 

The pad plane, with an effective area of 100 x 100 mm 2 , has 16 pad rows at a pitch of 6.3 mm, each 
consisting of 1.17 x 6 mm 2 rectangular pads arranged at a pitch of 1.27 mm. The neighboring pad 
rows are staggered by half a pad pitch. Pad signals are fed to charge sensitive preamplifiers located 
20 cm behind the bulkhead of the gas vessel. The amplified signals are sent to shaper amplifiers via 



twisted pair cables, and then processed by 12.5 MHz digitizers [27J . 

The chamber gases are mainly Ar-CF4-isobutane mixtures at atmospheric pressure and room 
temperature. The gas pressure and the ambient temperature (~ 20°C) are periodically monitored 
since they are not controlled actively. 

The prototype TPC is placed in the uniform field region of a superconducting solenoid without 
return yoke, having a bore diameter of 850 mm, an effective length of 1000 mm, and a maximum 
field strength of 1.2 T (see Fig. [5]). A pair of scintillation counters sandwiching the prototype TPC 
is used to trigger a data acquisition system upon detection of cosmic ray tracks. 

The measurement of gas gain was carried out using a dedicated small test chamber and an 55 Fe 
source without magnetic field. The small chamber has a GEM stack with a configuration similar 
to that of the prototype TPC and a 10-mm drift gap for photon conversion with a field strength of 
200 V/cm. Signal charge induced on its pad plane is read out by a combination of a preamplifier, a 
shaper amplifier and an analog to digital converter. 
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Figure 5: Photograph of the prototype ready to be inserted into the magnet. 



3. Gas gain measurements 

A small value of the diffusion constant (D) and a reasonably large effective number of electrons 
(iV e fj) are crucial properties of the LCTPC gas as mentioned in Section 1. Another major practical 
issue is the gas gain in the amplification gap. The demonstration of stable operation of GEMs at an 
adequate gain is important particularly in a gas containing CF4 since CF4 molecules dissociatively 
attach hot electrons during their passage through the GEM stack. We therefore first measured the 
gas gains in several Ar-CF4 based gas mixtures as function of the high voltage across each GEM foil 
(Vgem)- The electric fields in the transfer and the induction gaps were kept at 1.6 kV/cm throughout 
the measurements. 

The results are shown in Fig. [6]for Ar-CF 4 (3%) and Ar-CF 4 (3%) with 2% admixture of methane, 
ethane, propane or isobutanelll The gain obtained with a P10 gas (Ar-CH 4 (10%)) is also included as 
a reference. Figure [7] gives the measured gas gain as a function of the isobutane content. The signal 
could not be observed without isobutane at Vgem = 240 V. 

The gain study shows that the addition of a small amount of isobutane (or propane) to Ar-CF 4 
drastically lowers the operating GEM voltages which is required to obtain an adequate gas gain, as 
has already been found [13]. This is most likely due to the Penning effect since methane (ethane) 
is much less (moderately) effective as an additive gas@ In addition, Magboltz shows that a small 
content of isobutane reduces the electron attachment by CF 4 molecules in the transfer /induction 
gaps by cooling down the migrating electrons, thereby increasing the effective gain of a GEM stack 
(see Fig. [8]). It should be noted that high electric fields in the transfer and/or induction gaps could 
cause a significant loss of effective gain (even with isobutane). In the amplification region (E > 10 
kV/cm) the gas multiplication dominates. 

From the measurement above we decided to carry out the subsequent performance tests with the 
prototype TPC using a mixture of Ar-CF 4 (3%)-isobutane(2%), with Vgem = 240 V and the transfer 
and the induction fields of 1.6 kV/cm, resulting in an effective gas gain of ~ 8000. 



7 The CF4 concentration is hxed to 3% in the present study somewhat arbitrarily. It should be noted, however, 
that a higher CF4 content requires a higher drift field corresponding to the drift- velocity maximum. 

8 The first metastable excited state of argon is at 11.55 eV while the ionization potentials of methane, ethane, 
propane and isobutane are 12.75, 11.49, 11.07 and 10.78 eV, respectively. 
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Figure 6: Measured gas gain as a function of GEM Voltage. 
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Figure 7: Measured gas gain as a function of isobutane concentration. Gas: Ar-CF4(3%)-isobutane, Vqem = 240 V. 
The curve is only to guide the eye. 




Figure 8: Attachment rates (left) and effective ionization coefficient (right) as function of the electric field given by 
Magboltz. 
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Figure 9: Pad responses for different drift distances, increasing from the left to the right up to 257 mm. The upper 
(lower) plots are for B = T (B = 1 T). Gas: Ar-CF 4 (3%)-isobutane(2%), Drift field: 130 V/cm. See text for the 
definitions of the horizontal and vertical axes. The solid curves represent the gaussians fitted through the data points. 



4. Performance of prototype 

The data analysis of the recorded cosmic ray events was carried out using the traditional MultiFit 
(DoubleFit) program j28|. The hit coordinate on each pad row was determined by a simple charge 
centroid (barycenter) method using the pad signals within a charge cluster in the row, without any 
correction to the possible bias due to the finite pad pitch. Our major concern is the spatial resolution 
along the pad-row direction (erx) since the requirement to the resolution along the drift direction 
(er z ) is moderate for the LCTPC. The resolution o"z was measured to be 400 jttm - 700 /zm, depending 
on the drift distance and the drift-field strength. 

We present in this section the drift properties of the gas and the performance of the prototype 
TPC obtained using cosmic rays nearly perpendicular (± 2°) to the readout pad rows and parallel 
(± 8°) to the pad plane, unless otherwise stated. 

4-1. Drift properties 

Examples of the observed pad responses are plotted in Fig. 9, where the normalized charge on the 
pad (<&/ J2j Qji with g« being the signal charge on pad i ) is plotted against the distance of the pad 
center measured from the charge centroid, on a track by track basis. The data points are fitted with 
a gaussian (solid curve in the figure) in order to obtain the pad-response width for each region of the 
drift distance. Fig. 10 shows the pad-response width squared ( function of the drift distance 

(z) and a fitted straight line. The diffusion constants (D) were determined using the relation 

^PR = ^pro + D 2 -z . (3) 

The constant term is given by 

up 1 

^pro = 12 + ct p rf ' ( 4 ) 

where w denotes the pad pitch (1.27 mm) and oprf is the width of the avalanche-charge spread in the 
GEM stack along the pad row direction, for single drift electrons 0, [29|. The intersects of the fitted 
straight lines were always slightly larger than those calculated using Eq. (4) assuming the diffusion 
constants given by Magboltz in the transfer and induction gaps of the GEM stack to estimate oppj?- 
The slopes of the fitted lines give the diffusion constants, which are summarized in Table 1 along 
with the measured drift velocities. 
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The drift velocity was determined from the full width of the drift-time distribution and the 
maximum drift length (257 mm). Figure 11 shows an example of the drift time distribution. The 
errors of the measured drift velocities were estimated conservatively assuming reading errors of ± 1 
bin (80 nsec) both for the rising and falling edges of the distribution. 
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Figure 10: Pad-response width squared vs. drift distance. Gas: Ar-CF 4 (3%)-isobutane(2%), Drift field: 130 V/cm. 




Figure 11: Drift time distribution. Gas: Ar-CF 4 (3%)-isobutane(2%), Drift field: 100 V/cm, B=0T. 

The measured drift properties are plotted in Fig. 12 along with the corresponding Magboltz 
predictions, which closely reproduce the data points. The reliability of Magboltz for this gas mixture 
was thus confirmed for the practical range of the electric field for TPC operation, i.e. between those 
corresponding to the diffusion minimum and the drift-velocity maximum. 

4-2. Spatial resolution 

The spatial resolution along the pad row direction was obtained by plotting the geometric mean 
residuals against the drift distance on an event-by-event basis. A demonstration of the geometric 
mean method and the definition of the geometric mean residual are given in Appendix A. Figure 13 
shows examples of the scatter plots thus obtained and the width along the ordinate axis as a function 
of the drift distance. In order to estimate the widths (standard deviations) the scatter plot was binned 
by drift distance and fitted by gaussians. The figures clearly show the improvement of the resolution 
due to a 1 T axial magnetic field. It should be noted that the resolution obtained with the magnetic 
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Figure 12: Measured transverse diffusion constants and drift velocities. Gas: Ar-CF4(3%)-isobutane(2%). The 
Magboltz predictions (curves) are also shown for comparison. 



Table 1: Drift properties and N eS in Ar-CF 4 (3%)-isobutane(2%) 

(a) B = T 



E (V/cm) 


Drift velocity (cm//xs) 


Diffusion constant (/ihi/a/chi) 


(/im/ a/chi) 


N eS 


Measured 


Magboltz 


Measured 


Magboltz 


100 


4.45 ± 0.13 


4.50 


337 ± 12 


327.8 ± 2.4 


74 ± 2 


20.7 ± 1.8 


130 


5.67 ± 0.20 


5.67 


310 ± 11 


316.7 ± 1.9 


67 ± 2 


21.1 ± 1.9 


165 


6.61 ± 0.28 


6.69 


305 ± 8 


311.7 ± 2.5 


62 ± 2 


24.1 ± 1.5 


200 


7.38 ± 0.34 


7.35 


306 ± 8 


311.2 ± 2.2 


63 ± 2 


23.8 ± 1.4 


250 


7.74 ± 0.38 


7.81 


305 ± 10 


315.4 ± 2.4 


65 ± 2 


22.0 ± 1.8 


(b) B = 1 T 



E (V/cm) 


Drift velocity (cm//is) 


Diffusion constant (/ihi/a/chi) 


Measured 


Magboltz 


Measured 


Magboltz 


100 


4.52 ± 0.13 


4.50 


72 ± 2 


73.8 ± 0.7 


130 


5.74 ± 0.20 


5.67 


70 ± 2 


72.8 ± 0.6 


165 


6.69 ± 0.28 


6.69 


77 ± 2 


77.3 ± 0.6 


200 


7.30 ± 0.33 


7.35 


88 ± 2 


85.5 ± 0.7 


250 


7.65 ± 0.36 


7.81 


106 ± 2 


100.9 ± 0.7 
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Figure 13: Geometric mean residual vs. drift distance for B = T (left) and B = 1 T (right). Gas: Ar-CF4(3%)- 
isobutane(2%), Drift field: 130 V/cm. The widths of the residuals are plotted against the drift distance in the lower 
figures, with a fitted curve (Eq. (1)) for B = T and a calculated resolution for B = 1 T, respectively. 

field can not be fitted with Eq. (1) because of the degradation at small drift distances due to the 
finite pad-pitch effecfl The curve shown in the figure for B = 1 T is the result of the analytic 
calculation described in Appendix B and in Refs. (a Jiij]. 
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Figure 14: Spatial resolution vs. drift distance. Gas: Ar-CF4(3%)-isobutane(2%), Drift field: 250 V/cm. Data points 
for P5 gas [5J are also plotted for comparison. Gas: Ar-CH4(5%), Drift field: 100 V/cm. The curves in the figure are 
Eq. (1) fitted through the data points for B = T, and calculated resolutions for B = 1 T. 

In Fig. 14 the resolutions are compared to those measured previously with a P5 gas in the beam 
test [5| for drift fields corresponding approximately to the drift- velocity maxima. The figure shows 
the smaller degradation of the resolution at long drift distances due to diffusion in Ar-CF 4 -isobutane. 
The effective numbers of electrons (-/V e ff) f° r B = T listed in Table 1 were obtained readily from D 2 



9 This effect is much less prominent in the absence of magnetic field because of larger diffusion constants. 
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given by the pad response (Eq. (3)) and D 2 /N e s determined from the spatial resolution (Eq. (1)), as 
function of the drift distance. The average value of N e g for Ar-CF 4 (3%)-isobutane(2%) is 22.6 ± 1.5, 
which should be compared to 22 ± 2 obtained with a P5 gas [5]. All these values are comparable 
to an estimation It should be noted, however, that the apparent values of N e g could be smaller 
than the true values because of electronic noise and inadequate calibration of the readout electronics 
(see Appendix C). Therefore the measured values are the lower limits of the genuine N e g determined 
by the primary ionization statistics and the avalanche fluctuation. 

As mentioned above it is difficult to deduce the value of N ef[ using Eq. (1) for B — 1 T because 
of the finite pad-pitch effect unless the drift region is long enough. Actually the observed resolution 
is the quadratic sum of the diffusion contribution, the finite pad-pitch term and a constant offset 
(~ ffYn 1 10 l. The last term, which may depend on experimental conditions, can be estimated as the 
difference of the measured resolution squared and the calculated one. 

As an example, each component contributing to the resolution is shown in Fig. 15 along with the 
measured resolution for = 250 V/cm. In the calculation the diffusion constants in the drift region 
and in the transfer/induction gaps given by Magboltz were used. The value of N e ff was assumed 
to be 23. The constant offset is about (48 /im) 2 and is due most likely to electronic noise and 
inadequate calibration of the readout electronics. Possible major sources of the constant term are 
briefly discussed in Appendix C. 
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Figure 15: Spatial resolution squared vs. drift distance. Gas: Ar-CF4(3%)-isobutane(2%), Drift field: 250 V/cm, B 
= 1 T. The curves are the calculated resolution squared (solid line) and its components: the finite pad-pitch term 
(dashed line) and the diffusion term (dash-dotted line). The horizontal straight line represents the fitted constant 
offset. 

The net diffusion contribution can thus be extracted by subtractions and is shown in filled squares 
in the figure. By fitting a straight line through those data points, N e g under the 1 T magnetic field 
could be estimated. The resultant values are about 21 and 23, depending on the diffusion constant 
assumed, the Magboltz value (100.9 /xm/ 'y/cm) or the measured value (106 /im/y 7 ^^. The values 
of iV e ff thus obtained is consistent with those found from the resolutions in the absence of magnetic 



10 Strictly speaking, the constant term (<7xo) is slightly greater than the constant offset since the diffusion term 
contains small but finite contribution to the constant term (see Eq. (|B.5|) in Appendix B). 

11 The diffusion constant and the value of N e g assumed in the calculation do not affect the resultant N e g since they 
are used only to estimate the constant offset. 
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field (see Table 1), indicating that the effective number of electrons is not affected by the existence 
of the magnetic field as well as by the variation of the drift field. 

5. Expected resolution of an LCTPC 

Our measurements of the drift properties in Ar-CF 4 (3%)-isobutane(2%) are well reproduced by 
the Magboltz code for B = and 1 T. This indicates that the cross sections of the gas molecules 
involved are correctly incorporated in the code at least for the corresponding electron energies for 
Ar-CF 4 -isobutane seen in Fig. 3 and that the influence of the magnetic field on the transverse 
diffusion is properly included. Since the energy (random velocity) distribution of drifting electrons 
is not affected by the existence of axial magnetic field Magboltz is expected to give correct values of 
diffusion constants under stronger magnetic fields. In addition, the effective number of electrons is 
measured to be about 23 for the pad row pitch of 6.3 mm, which is affected little by the magnetic field 
if any. Therefore it is now possible to estimate the spatial resolution of the LCTPC using Eq. ( IB. 41 ) 
in Appendix B. 

Figure 16 shows the expected resolutions as function of the drift distance up to 2.5 m for a GEM- 
equipped LCTPC. In the calculation of the resolutions the transfer gaps and the induction gap of 
4 mm are assumed for a triple GEM stack, with the electric fields in the gaps of 1600 V/cm while 
the drift field is set to 250 V/cm corresponding approximately to the drift- velocity maximum. The 
resolution better than 100 /im is thus feasible throughout the sensitive volume for tracks perpendicular 
to the pad row under an axial magnetic field of 4 T, at which the diffusion constant is predicted to 
be 29 /im/ \/cm by Magboltz, if the constant offset is kept reasonably small. 




Drift distance (m) 



Figure 16: Expected spatial resolution of a GEM-equipped LCTPC for tracks perpendicular to the pad rows. Gas: 
Ar-CF4(3%)-isobutane(2%), Drift field: 250 V/cm, B = 4 T. The transfer and induction gaps are assumed to be 4 
mm. 

6. Conclusions 

We have measured the diffusion constant, the drift velocity and the spatial resolution, as well 
as the gas gain, using a GEM-equipped TPC prototype operated in a gas mixture of Ar-CF 4 (3%)- 
isobutane(2%) under magnetic fields of T and 1 T. The prototype has been working stably with a 
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sufficient gas gain in the cosmic ray test and gave a reasonable effective number of electrons (N e g), 
indicating a small loss, if any, of drift electrons at the entrance to the GEM stack. A reasonable 
value of N ef[ means also that the gas-gain (avalanche) fluctuation is moderate. A small amount of 
isobutane (~ 2%) was found to be very effective to lower the operating high voltages of the GEMs, 
most likely due to the Penning effect. 

The measured diffusion constants (D) in Ar-CF4(3%)-isobutane(2%) are consistent with those 
given by Magboltz both for B = T and B — 1 T, demonstrating the reliability of the Magboltz 
predictions for gas mixtures containing argon, CF 4 and isobutane. This gas mixture is confirmed to 
be cooler and faster than conventional TPC gases such as P5 and to ensure better azimuthal spatial 
resolution at long drift distances under a high axial magnetic field. 

The spatial resolution better than 100 /xm, the target of the LCTPC, is expected to be within 
reach with a GEM-based TPC operated in Ar-CF 4 (3%) plus a small amount of isobutane under a 4 
T magnetic field. 

Appendix A. Geometric mean method 

In this appendix we give a simple demonstration of the geometric mean method applied in the 
analysis to estimate the spatial resolution. 

First, in the case where the hit point in question (say, Xi, on the i-th pad row) is excluded in the 
track fitting, the residual is given by 



where Xi represents the estimator for the i-th point given by the track fitting using the remaining 
hit points. Since X{ and X{ are statistically independent the variance of the residual is given by 



the sum of the true spatial resolution and the tracking error. 

Next, in the case where the hit point in question is included in the track fitting, the estimator for 
the i-th hit point is the weighted mean of xi and xf. 



with w Xi (w Xi ) being the corresponding weight: 1 / cr? i (l/cr^.). The residual is hence given by 





= a r . + (Ta 



(A.l) 



Ax[ = Xi - x\ 




{Xi Xi) 



■ Axi . 



(A.2) 



The variance of the residual in this case is therefore 
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> 2 } 



cxcl 



(A.3) 



Combining Eq. flA.ip and Eq. ( 1A.3jl . the true spatial resolution is given by 



0~ Xi = V^cxcl ■ fincl • (A.4) 



Equation (IA.4j) demonstrates that the spatial resolution is obtained by taking the geometric 



mean of the two widths: one for the Axi distribution and the other for the Ax- distribution. 
One may instead plot, on an event-by-event basis, geometric mean residuals defined as Ax" = 



S ■ J (xi — Xi) ■ (xi — x'i) , with S being the sign of x% — Xi (or equivalently Xi — x[ from Eq. ( IA.2j) ) 
and interpret its width as the spatial resolution since 



a Ax'/ — \S 2 ' (Xi — Xi) ■ (Xi — Xj) 

- (s ■ \J (xi - Xi) ■ {xi - x~^)\ 

{((Xi - Xi) 2 ") - (Xi - x^ 2 } 



Appendix B. Analytic expression of the spatial resolution 

One way to estimate the spatial resolution of a TPC along the pad-row direction is to write a 
realistic Monte-Carlo simulation code. This technique is applicable to any complicated situation, and 
has been developed by several groups. On the other hand, an analytic approach is applicable only to 
a restricted case where incident particles are normal to the pad row. However, the resultant formula 
is rather simple and is sometimes enlightening as shown below. Though a numerical calculation is 
needed to evaluate the formula, the demanded CPU time is much less than a Monte-Carlo simulation. 
In addition, the analytic calculation can be used to check the reliability of a Monte-Carlo simulation 
program, which is usually long and complicated. This appendix is devoted to briefly summarize our 
analytic approach, based on the following assumptions. 

1. Particle tracks are normal to the pad row. 

2. Track coordinate is determined by the charge centroid (barycenter) method. 

3. Displacement of arriving drift electrons due to E x B effect near the entrance to the detection 
device is negligible. 

4. Displacement of arriving electrons due to the finite granularity of amplification elements of the 
detection device (line intervals in MicroMEGAS or a hole pitch in GEM) is negligible. 

Appendix B.l. Electronic noise contribution 

Let us first consider the electronic noise contribution to the resolution degradation. The charge 
centroid of the signal including the noise is given by 
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where qj (g'- ) is the signal (noise) charge on pad j and x* is the coordinate of the center of pad 
j. For simplicity we assume that J2jQj ^ J2jQj, therefore J2j(Qj + Qj) ~ J2jQj, an d that (<^) = 0, 
(Qj ' Qk) = (q' 2 ) ' djk an d (Qj ' Qk) = 0- Then, with x being the true track coordinate, the squared 
resolution along the pad-row direction is given by 



<4 = ((x-xY) 
1 



Eifaj + Qjl 



. 2 

+ g')(a;* -x 



1 j 
1 ^ 



Hfe + Qj) (xj - x) + + fyXg* + Q k ){Xj - x)(x k - x 

3 \H 3 -lj I ' ^ j j/k 

J2($ + Q?)( x *j-x) 2 + I2Wk{x*j-x)(xl-x] 

J v .? j/k 

1 \ 2 \ // 1 \ 2 



E,- <?; 



Ej Qj j J / \\ Yjj Qj 

The second term represents the noise contribution whereas the first term gives the resolution obtained 
in the ideal case in which the electronic noise is absent. 

Let us define Qi as the total gas-amplified signal charge on the pad row created by drift electron 
i. The total charge Qi is given by J2j Qji, where is the signal charge contribution of drift electron 
i, on pad j (qj = J2i Qji)- Then J2j Qj — Ej Ei Qji — Ej Qi- Assuming that the total number of drift 
electrons per pad row (JV, not a constant) is large enough 

N 

Y,Qj = Y,Qi~ N ■ (Q) > 

j i=l 

where (Q) is the average total signal charge per drift electron. 
Now the noise contribution is given by 

4) noi se-(^)-^-E(^-^) 2 - (B.2) 

This term increases with the signal charge spread along the pad row direction given by Eq. (3) 
because of the increase of the active number of pads (the summation range of j around x). It should 
be noted that be written as j ■ w with w being the pad pitch. 

Appendix B.2. Analytic evaluation of spatial resolution 

In this section we evaluate the first term in Eq. (IB.lj) . Let us redefine the charge centroid as 
follows. 

_ Ej Ej=l Qji ' x j _ Ej=l Qi Ej Fj( x i) ' Xj 
Ej Ei=l Qji Ei=l Qi 

where (see Fig. IB. If) 

N : total number of drift electrons per pad row, 
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X.: 



arrival position of electron i at the entrance to the detection device, 



e j ■ w : central coordinate of pad j (j = ■ 

J jw—w/2 

/(£) : (normalized) pad response function (PRF). 



•,-2,-l,0,+l,+2, 



It should be noted that the PRF here is determined solely by the charge spread in the detection gap 








Pad; 






Pad response 




function 





Readout pads 



Figure B.l: Illustration of the relevant variables. 



and does not depend on the pad pitch whereas conventional PRFs do. 
Ther£! 



a 2 x {x) = ((X-x) 2 ) 



i i - Xj 



J2i=i Qi 



— X 



In what follows, for an arbitrary function /, 





3 (n/^( 






(HQ)) 





N 



N 



where P x {%i) and Pq(Qi) denote the probability density functions, respectively for the arrival position and the signal 
charge of drift electron i whereas N is the total number of drift electrons per pad row (not a constant, actually). It 
should be noted that (Ii(x) ■ h{Q)) — (h(%)) x ■ (h(Q))n since X{ and Qi are uncorrelated. Note also that the arrival 
positions of drift electrons (xj) are not correlated. 
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Ej=l Qi Ej ' ( x j x ) 



(E» Qt) : 



(E» Q0 : 



eo?-(e^o-(^-*)) 2 

i 3 

+ E Q<Q« • (E ^te) • (4 - 5) E F M) -(x* k -x 

ijtl ^ k k 

(e^(e^-^-z) 2 

+ E QiQi ■ ( (E • 4 - ^) (E 
(EQ?-((E^)^-s) 



X J., X 



+ (E QiQi - E Q') ■ ( E ^ (*) • 4 - ^ 



\ - / v n 2 



E • - ^ + (^) g • ((E W • • ■ *t) x - (E F i( x) 



X 3 

3,k " ' 



E<^)) • ^ - 5) 2 + ((§^5) • (E<**(*) • • 3 • 4 - (E(*K*)> • xf 2 
E<*K*)> • - 5) 2 + 1 • |p • (E W*) • F ^)> • 3 • 4 - (EW*)> • ^*) 2 ) 



Averaging over N and substituting j ■ w and k ■ w, respectively for x* and x* k we get 
o*x@) « (E> • " + ^ • fEi^ 2 • <*i(*) • - (E> • (^)>V ) , (B.3) 



where 



1 _ / 1 \ (Q 2 ) 1 1 \ A . 4 



TT '^4= TT ' 1 



2 / ' 



iVefr W/ (q> 2 \iv/ V (qy 

/oo 
P x (x) • Fj(x) dx , 
-oo 

/oo 
P x (x) ■ Fjix) ■ F k (x) dx , 
-oo 



with 

. . 1 / (x — x) 2 

pjw+w/2 

Fj(x)= f(S-x)dt, 

J jw—w/2 

The parameter cr<f in the definition of P x (x) denotes the (lateral) standard deviation of the arrival 
position of drift electrons at the entrance to the amplification device measured from x, and is given 
by D ■ y/z. 
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Hence the spatial resolution as a function of the drift distance (z) can be numerically evaluated 
once the pad pitch (w), the diffusion constant (D), the pad response function (/), and the effective 
number of electrons (N eS ) are given. In our case the pad response function is assumed to be a 
gaussian 

/(£) = _L-exp( L 



with a a being the charge spread in the GEM stack. 

Though Eq. flB.3j) may appear a little complicated its meaning is quite simple: the first term 
is the bias inherent in the charge centroid method (finite pad-pitch term) while the second term 
represents the variance around the bias, divided by N e $ (diffusion term). This is shown clearly by 
rewriting Eq. (1B.3j) as 

4(x) « feFjix) -x*-x)\ ■ x* - feFfa) ■ x] 

It should be pointed out here that a\ depends on the position of x relative to the corresponding 
pad center, and that the beam spot size is usually much larger than the pad pitch. Therefore unless 
the incident positions of incoming particles are measured precisely by external trackers (e.g. by a 
set of silicon strip detectors) on an event-by-event basis, u x (x) obtained above (Eq. (1B.3j) ) has to be 
averaged over x in a range, say, [-w/2, -Ho/2]: 

1 r+ w / 2 
cr x = — ■ / &x(x) dx . 

W J -w/2 

In order to see the asymptotic behavior of the resolution it is convenient to rewrite Eq. flB.3j) : 



2 „ . x * _ ^ + _L ^ Fj{x) . x * -if)- (E F,(X) ■ X* - x) 2 



a x [x 



and define 



where 



1 r w / 2 1 
o\ = -- a 2 x (x) di ^ B 2 + -- ■ (A 2 - B 2 ) , (B.4) 



W J -w/2 AV 7 N, 



cS 



1 r w l 2 I \ 2 

1 rw/2 If \2\ 

* - --/ ((e *;-*))««*■ 



The term B 2 vanishes in the limit cr^ — > oo (z — > oo) since J2j (Fj{x)) x* — > x for any given x 
while the term A 2 can easily be shown to approach B 2 + o\ = B 2 + D 2 ■ z in the same limit, where 
<r d =o = B 2 } z=0 . 



£? 2 = B 2 L J=n = 5 2 L =n . Therefore 



4 -> — - • (B 2 + D 2 -z) (B.5) 

at long drift distances. 

Figure IB. 21 shows an example of the calculated resolution in our case. In the calculation the 
following values were assumed: w = 1.27 mm, o a = 279 /iin, D = 100.9 fim/y/cm (Magboltz value), 
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and N e g = 23. a a was estimated using the diffusion constant in the transfer and the induction gaps 
given by Magboltz. The dashed straight line shows the asymptotic behavior at long drift distances 
(Eq. flB.5j) ). which has an offset of Bq/N c r « (23 yum) 2 . This is the fundamental lower limit of the 
constant contribution to the spatial resolution (cr xo ), determined by the pad pitch (w), the charge 
spread in the detection gap (cr a ), and the effective number of electrons (N e g). 




Figure B.2: Calculated resolution squared for Ed = 250 V/cm and B = 1 T. See text for the dashed line. 



Appendix C. Contributors to the constant term 

In this appendix we qualitatively discuss the origin of the constant term contained in the spatial 
resolution (cr xo ) in our case. Possible main sources of the constant term are listed below for tracks 
perpendicular to the pad-row direction. 

1. Intrinsic track width (5-rays) (3o| . 

2. Multiple Coulomb scattering. 

3. Mechanical imprecision. 

4. Influence of the finite hole pitch of the GEM foils. 

5. Influence of the finite pad pitch. 

6. Electronic noise. 

7. Inadequate calibration of the readout electronics, such as channel-to-channel (pad-to-pad) in- 
equalities of the gain and inappropriate pedestal subtractions from the raw signals. 

The sources (1) and (2) are intrinsic and unavoidable while the others can be made sufficiently small 
in principle. 

The intrinsic track width certainly contributes to the constant term. The value of <t xo obtained 
with the magnetic field (~ (50 /im) 2 ) is significantly smaller than those for the resolutions obtained 
without the magnetic field (~ (100 yum) 2 , see Figs. 13 and 14 for example) because of the suppression 
of the transverse intrinsic track width by the axial magnetic field. The magnetic field reduces the 
effective range of 5-rays along the pad-row direction. 
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The mechanical precision such as the readout pad arrangement is expected to be good enough. 
The influence of the hole pitch of the GEM foils (horizontal pitch of 140 jum) was also found to be 
small by means of a Monte-Carlo simulation (see, for example, Fig. 4 of Ref. j^j). The constant 
term (at long drift distances) is caused by the finite pad pitch as well (see Appendix B, Eq. ( IB. 51) ). 
However, its contribution is small in our case, (18 fim) 2 - (23 /im) 2 , depending on the magnetic field 
(0 T - 1 T). 

The contribution of the angular pad effect for the azimuthal angle cut of ± 3° was estimated to 
be ~ (20 /im) 2 by a Monte Carlo simulation. Therefore the main contributors to a\ Q are expected 
to be (6) and (7), except for the unavoidable contribution from (1) and (2). 

The sources (6) and (7) have contribution increasing with the drift distance (z) as well as a 
constant term at z = because of the increase of the number of active pads used to determine the 
charge centroid, due to diffusion (see Figs. 9 and 10, and Eq. (IB.2j) for the contribution of electronic 
noise). Therefore they can affect the apparent value of iV e fj as well. 
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